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Ah st  me  x 

This  report  considers  an  array  of  simultaneous  under¬ 
water  nuclear  explosions  generating  piano  tidal  waves  on 
the  continental  shelf.  Computer  simulations  are  used  to 
study  this  wave  generation  process.  They  are  based  on  non¬ 
linear  shallow  water  theory  with  an  initial  condition  of  a 
stationary  raised  cone  of  water.  The  sensitivity  of  the 
wave  height  to  different  yields,  combination  of  different 
yields,  spacing,  alignment,  timing,  and  bottom  slope  is 
determined.  The  reference  values  used  for  this  study  are  a 
yield  of  five  kilotons  TNT,  one  kilometer  spacing,  water 
depth  of  100  meters,  and  no  bottom  slope.  In  addition  to 
the  sensitivities,  the  amount  of  non-linear  effect  is  illus¬ 
trated  by  comparison  of  wave  contour  printer-plots  generated 
by  linear  and  non-linear  shallow  water  theory.  The  computer 
results  indicate  the  wave  generation  is  relatively  insensitive 
to  realistic  variations  in  the  parameters. 
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PI.ANL  TIDAL  KAVF.S  CLNIRA'l  Li)  BY 
AN  ARRAY  OF 

SIMULTANEOUS  UNDERWATER  EXPLOSIONS 

I .  T  n  t  roduc t i on 


Backg  round 

This  thesis  originated  from  the  A  FIT  CNL-80M  class 
design  project  which  analyzed  a  shallow  underwater  mobile 
(SUM)  basing  alternative  to  the  land  based  MX  intercontinental 
ballistic  missile  weanon  system  (Ref  1).  The  SUM  system  would 
be  deployed  over  the  continental  shelf  of  the  United  States. 

This  large  area  of  deployment  with  the  location  of  the  SUMs 
unknown  would  require  a  large  number  of  warheads  to  defeat  the 
weapon  system.  However,  an  easy  defeat  of  the  SUM  based  weapon 
system  lias  been  proposed  which  uses  a  plane  tidal  wave  generated 
by  an  array  of  simultaneous  underwater  explosions.  The  tidal 
wave  would  travel  across  the  continental  shelf,  disabling  the 
SUMs.  Computer  simulations  based  on  shallow  water  theory  verify 
the  formation  of  a  plane  tidal  wave.  The  wave  energy  resulting 
from  the  underwater  explosion  was  approximated  with  an  initial 
condition  of  a  stationary  raised  cone  of  water.  Further  study 
of  this  plane  tidal  wave  generation  is  the  topic  of  this  thesis. 

Problem  and  Scope 

The  primary  purpose  of  this  study  is  to  determine  the 
sensitivity  of  the  plane  tidal  wave  generation  to  the  explosions' 
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yield,  yield  combinations,  spacing  between  explosions,  align¬ 
ment  of  explosions,  timing  of  explosions,  and  slope  of  the 
continental  shelf.  The  reference  values  used  are  five  kilotons 
TNT  for  yield,  one  kilometer  spacing  between  explosions,  and 
a  water  depth  of  100  meters.  The  yield  is  varied  from  2  to  8 
kilotons.  The  combination  of  two  yields  is  varied  as  a  ratio 
of  one  yield  to  the  other  yield  from  .25  to  1  with  an  average 
yield  of  5  kilotons.  The  spacing  is  varied  ±40oo  from  the 
reference  value  of  one  kilometer.  The  misalignment  is  varied 
from  0.0  to  0.5  kilometer  and  the  timing  is  varied  from  0.0  to 
5  seconds.  The  bottom  slope  is  varied  from  0.0  to  0.002. 

In  addition  to  the  sensitivities,  the  amount  t ha t  the  non¬ 
linear  effects  contribute  to  the  generation  of  the  plane  tidal 
wave  is  illustrated. 

Assumpt  ions 

There  are  two  main  assumptions.  The  first  is  that  an 
explosion  on  the  ocean  floor  which  produces  a  large  enough  gas 
bubble  for  blowout  will  generate  shallow  water  waves.  The 
second  assumption  is  part  of  shallow  water  theory.  It  assumes 
the  water's  vertical  acceleration  doesn't  affect  pressure  and 
t he  water's  vertical  velocity  doesn't  affect  the  conservation 
of  momentum  equation. 

Approach  and  Sequence 

The  general  approach  starts  by  analyzing  the  problem 


physically  to  determine  the  type  of  wave  generated 


real i s  t i c 


parameter  reference  values,  and  modeling  of  initial  condition. 
Next,  the  differencing  scheme  is  developed  by  deriving  the 
shallow  water  equations  and  differencing  t lie  equations.  The 
scheme  is  used  to  simulate  the  wave  generation.  The  wave 
heights  are  evaluated  to  determine  sensitivity  to  several  of 
the  parameters.  Finally,  the  non-linear  effects  are  illus¬ 
trated  . 


I  [  .  Analysis  o  1  1’  rob  1  cm 

The  analysis  of  this  problem  focuses  on  three  areas. 

First  is  the  rationale  and  basis  for  using  shallow  water 
theory.  Next  is  the  determination  of  the  reference  values 
for  the  independent  variables  and  their  range  of  variation. 

The  last  area  is  the  initial  condition  and  its  relationship 
to  the  explosions'  yield. 

Rat  ion- le  for  Shallow  Water  The or v 

The  rationale  for  using  shallow  water  theory  is  primarily 
based  on  the  type  of  cavity  formed.  Shallow  water  theory  can 
he  used  if  the  total  depth  of  the  water  moves  horizontally 
together.  The  explosions  used  for  this  problem  generate  a 
gas  bubble  witli  a  radius  greater  than  the  depth  of  water. 

Since  the  expanding  bubble  will  horizontally  excite  the  total 
depth  of  the  water,  shallow  water  waves  are  expected.  On  this 
basis,  shallow  water  theory  is  used. 

Research  at  l.os  Alamos  Scientific  Laboratory  (LAST)  found 
that  shallow  water  theory  was  not  applicable  for  a  near  surface 
explosion  (Ref  b) .  However,  the  difference  in  burst  location 
can  account  for  this.  The  near  surface  burst  used  by  LASL 
formed  a  hemispherical  surface  cavity  with  a  radius  equal  to 
1/6  the  depth  of  the  water.  The  surface  cavity  would  not  hori¬ 
zontally  excite  the  total  depth  of  the  water  at  the  same  time. 
Therefore,  shallow  water  waves  are  not  expected,  which  is 
consistent  with  shallow  water  theory  not  accurately  duplicating 
the  experimental  data. 
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Parameters  Re f e rcnc e  Values 

The  reference  values  of  1  lie  independent  parameters  and 
range  of  variation  were  selected  based  on  many  factors.  The 
water  depth  and  bottom  slope  are  based  on  the  average  depth 
and  slope  of  the  continental  shelf.  The  yield  was  selected 
based  on  wave  generation  efficiency  and  the  minimum  yield 
to  ensure  blowout.  The  spacing  was  limited  by  the  accuracy 
of  the  bomb  placement.  The  alignment  and  timing  are  assumed 
perfect  for  the  reference  value. 

Depth  and  Slope.  The  reference  water  depth  is  about  the 
average  of  the  continental  shelf's  depth  (  0  -  000  ft).  '1  he 
water  depth  reference  value  is  100  meters  with  no  variation 
except  when  bottom  slope  is  used.  The  average  bottom  slope 
of  the  continental  shelf  is  .0002  to  .0004  based  on  a  width 
of  300  to  600  miles.  The  reference  value  for  bottom  slope  is 
zero  with  a  variation  of  0.0  to  .002. 

liming  and  Alignment.  The  reference  value  for  timing  is 
zero  seconds  with  a  range  of  variation  of  0  to  S  seconds. 

The  reference  value  for  alignment  is  also  zero  with  a  range 
of  variation  of  0  to  500  meters.  Misalignment  distance  is 
the  perpendicular  distance  from  the  array  centerline  to  the 
misaligned  explosion. 

Spacing.  The  spacing  reference  value  is  limited  at  the 
lower  end  by  the  accuracy  of  the  warhead  placement.  According 
to  Scientific  Amer i can ,  the  1085  estimated  circular  error 


probable  (('HIM  lor  "ussian  missiles  is  ..12  kilometer  (Ref  2  :  54 )  . 
The  reference  spacing  was  limited  to  a  distance  of  4*CHP. 

Based  on  this,  the  reference  value  of  one  kilometer  is  used 
with  a  variation  of  .  (>  to  1.1  kilometer. 

Yield.  The  reference  yield  (Y)  is  based  on  the  efficiency 
of  wave  formation  and  the  minimum  yield  required  for  a  blowout 
case.  The  efficiency  as  a  function  of  yield  is  indicated  by 
the  wave  height.  The  wave  height  (h)  is  proportional  to 
(Ref  5:133).  Therefore,  the  lowest  possible  yield  will  be  used 
to  obtain  the  highest  efficiency. 

The  minimum  yield  required  to  produce  a  maximum  hemisper- 
ical  bubble  with  a  radius  equal  to  t he  water  depth  of  100  meters 
is  1 . 7  kilotons.  This  is  based  on  the  following  approximations 
and  calculations.  The  energy  required  to  form  the  hemispherical 
bubble  is  approximated  by  the  energy  required  to  form  an  equal 
volume  spherical  bubble  with  coinciding  centroids.  Approximately 
half  of  the  explosion's  energy  is  absorbed  by  the  ocean  floor 
(Ref  5:12b).  Using  these  approximations  and  hq  (1)  for  the  maxi¬ 
mum  bubble  radius  (Ref  5:113),  the  minimum  yield  is  determined 
as  follows.  1,  is  the  max  bubble  radius,  Y  is  the  yield,  1) 
is  the  depth  of  burst,  R  is  the  radius  of  the  spherical 
bubble,  and  R^  is  the  radius  of  t he  hemispherical  bubble. 
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To  obtain  the  minimum  yield,  use  Hq  (3)  with  I.  =  R .  , 

D  equal  to  the  depth  of  the  hemisphere's  centroid,  and 
Y  =  i  Y  to  account  for  the  energy  absorbed  by  the  bottom. 

j  Y ( kt )  =  (62.5  +  10) 

Y  =  1.7  ktons 


Tlierefore,  the  yield  should  be  equal  to  or  greater  than  the 
minimum  yield  of  1.7  kilotons.  A  reference  yield  of  5  kilotons 
was  selected  since  it  is  small  and  vet  more  than  twice  the 
minimum  yield  to  ensure  shallow  water  waves. 

Initial  Condition 

The  initial  condition  used  to  represent  the  wave  energy 
imparted  by  the  explosion  is  very  simplistic.  It  consists  of 
a  raised  stationary  cone  of  water.  The  cone  surface  has  a 


constant  slope  of  0.2  with  its  volume  dependent  on  the 
yield. 

The  relationship  between  the  yield  and  the  potential 

energy  of  the  cone  is  based  on  the  scaling  laws  for  the 

blowout  case  (Ref  5: 132-1213')  .  In  the  blowout  case,  the 

wave  height  is  scaled  by  Y^*'  .  Since  the  potential  energy 

of  the  wave  is  proportional  to  the  wave  height  squared,  it 

is  scaled  by  Y^3  .  Therefore,  the  potential  energy  (PH) 

1/3 

of  the  cone  is  equal  to  some  constant  times  Y 

PH  =  Const*Yl/3  (4) 

c  on  e 

The  constant  is  obtained  by  analyzing  the  potential 

energy  of  the  cavity.  Starting  with  the  minimum  yield  of 

1.7  kilotons,  the  maximum  bubble  radius  (1.)  is  100  meters. 

Assume  the  water  displaced  by  the  bubble  is  now  located  at 

the  surface.  The  amount  of  increase  in  potential  energy  is 

determined  as  follows:  where  p  is  water  density  (1000  kg/m  ) 

2 

and  g  is  the  acceleration  of  gravity  (9.8  m/s*'}  . 

APT.  =  (62.5)(j  £  n  L3)  pg  =  1.28*1012  Joules  (5) 

As  the  bubble  collapses,  much  of  the  potential  energy  is 
lost  in  turbulence.  One  half  of  the  wave  energy  is  assumed 
to  be  lost  in  turbulence.  Based  on  this  assumption,  the  constant 
is  determined  as  follows: 
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PE  (J)  =  5 . 4  *  1 0 1 1  Y1/3(kt) 

cone 
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The  size  of  the  cone  is  determined  by  the  amount  of  wave 
energy  imparted  by  the  explosion.  The  potential  energy  of 
the  cone  of  water  is  set  equal  to  the  energy  obtained  from 
Eq  (7).  Solving  for  the  height  (h)  of  the  cone's  peak 
determines  the  cone's  size  with  a  constant  slope  (s)  of  0.2. 
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III. 


(.1  output  or _ Prop,  ram  Dove  1  opine nt 

The  computer  propram  development  includes  the  derivations 
and  differencing  of  the  equations  and  treatment  of  t lie  boundary 
conditions.  Both  the  non-linear  and  linear  equations  will  be 
derived  and  differenced.  fn  addition,  a  check  of  the  differ¬ 
encing  scheme  is  accomplished. 

Derivation  of  r.quations 

The  shallow  water  wave  equations  are  derived  from  the 
conservation  of  mass  and  conservation  of  momentum  equations 
using  the  3 -dimens iona 1  rectangular  coordinate  system.  The 
viscosity  and  compressibility  of  the  water  are  ignored.  The 
coordinate  system  is  as  shown  in  Figure  1  with  the  free-standing 
water  surface  at  z  =  0  .  The  actual  water  surface  position 

is  described  by  H(x,y,t)  and  the  bottom  is  described  by 
b(x,y)  .  The  averaged  velocity  components  u  ,  v  ,  and  w 
are  in  the  x  ,  y  ,  and  z  directions,  respectively,  with  S 
representing  the  velocity  vector  (u,v,w) 

The  derivation  starts  with  the  conservation  of  mass  equation 
(Ref  4:2,24)  and  its  integration  with  respect  to  z 
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Now  the  first  approximation  must  be  made  for  the  shallow 
water  theory.  The  vertical  acceleration  of  the  water  is 
assumed  not  to  affect  the  pressure  (P)  (Ref  8:23-24).  This 
results  in  hydrostatic  pressure,  I*  =  up(ll-z)  where 
h  -  z  ‘l  H  .  Hydrostatic  pressure  provides  a  horizontal 
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pressure  gradient  which  is  i  tulependen?  of  z  .  Therefore, 
the  horizontal  accelerations  and  velocities  are  independent 


of  z  .  With  this  approximation,  the  integration  of  the 
conservation  of  mass  equation  (17)  is  completed. 
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The  conservation  of  momentum  equation  (Kef  4:3,24')  is 
used  to  derive  the  velocity  equations.  The  previous  approxi¬ 
mation  is  expanded  to  accomplish  this.  The  vertical  velocity 
w  is  assumed  to  have  a  negligible  effect  on  the  conservation 
of  momentum.  Therefore,  the  velocity  vector  S  becomes  (u,v)  . 
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The  equation  for  v  is  obtained  similarly. 


? 

av  tv'  ^  auv  .  an 

at  av  ax  K  av 


o 


(20) 


(21) 


(22) 


1  3 


Equations  (19) ,  (2L)  and  (221  are  the  non-linear 
shallow  water  equations.  The  linear  shallow  water  equations 
are  obtained  by  assuming  u  ,  v  ,  I!  ,  and  their  derivatives 
are  small  enough  so  their  squares  and  products  can  be  neg¬ 
lected  (Ref  8:24).  This  results  in  linear  equations,  where 
-b  is  simply  the  depth  of  the  water. 


an 

3 1 


+ 


3 

3  x 


[u( -bj 


+ 


_3_ 

3y 


Iv(-b)  | 


3u 

3 1 


3  v 

at 


+ 


+ 


o 


0 


0 


(23) 


(24) 


(25) 


If  the  water  depth  is  assumed  constant,  the  variables 
u  and  v  can  be  eliminated  as  follows  to  obtain  the  linear 
wave  hq  (30)  (Ref  8:25). 
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This  shows  that  the  characteristic  wave  velocity  is  cg(-bj 


Difference  Scheme 


The  second  order  explicit  difference  scheme  is  three 
levels  in  time  and  space,  conservative,  and  has  an  amplifi¬ 
cation  factor  of  unity.  Correct  space  and  time  centering 
provide  an  amplification  factor  of  unity  if  the  Courant- 
Fr iedr ich-Lewy  criterion  is  satisfied  (Ref  7:288).  This 
criterion  is  |u||-|  +  |v|~|  1  1  .  Using  the  charac¬ 

teristic  velocity  equally  distributed  between  u  and  v  , 
the  criterion  value  is  0.44  for  At  =  .5  sec  and  Ax  =  Ay 
=  50  meters  .  The  difference  equations  are  based  on  the 
integral  conservation  of  mass  lit}  (19)  and  the  conservation 
of  momentum  Hqs  (21)  and  (22). 
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Kinematic  viscosity  can  be  added  to 
equations  to  account  for  viscosity  or  to 
instability.  The  kinematic  viscosity  is 


the  last  two 
help  prevent 
incorporated  by 
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adding  to  the  right  side  of  I'.qs  (.>_’)  and  (ad)  a  kinematic 
viscosity  coe f f  ic  i  en t  ,  multiplied  h\  the  Laplacian  of  u 
and  v  ,  respectively.  A  slight  diffusion  coupling  might 
also  be  added  to  keep  the  independent  conservative  flows 
in  step.  This  is  done  hv  adding  the  product  of  a  coupling 
coefficient  and  the  Laplacian  of  II  to  the  right  side  of 

l:q  (31). 

The  difference  equations  use  the  indices  i  ,  j  , 
and  n  for  the  variables  \  ,  y  ,  and  t  ,  respectively. 
The  equations  are  second  order  in  time  and  space  with  no 
kinematic  viscosity  or  diffusion  coupling  as  shown  below. 
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The  linear  difference  equation  is  based  on  the  linear 
wave  l!q  (30).  It  is  second  order  with  space  and  time 
centering. 
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The  difference  scheme  used  a  standard  rectangular  mesh 
with  all  the  point  values  calculated  for  each  time  step. 
Symmetry  was  used  as  much  as  possible  with  the  aid  of 


reflecting  boundaries. 


Bound  a  r  y  Cuiul  i  l  1  on  s 

Reflecting  boundary  com!  i  t  i on s  arc  used  in  tlic  program 
to  take  advantage  of  symmetry  (Ref  7  : 2 8 X  -  J 0  1  )  .  With 
reflecting  boundaries,  the  boundary  conditions  are  set  up 
to  ensure  no  mass  flow  across  the  boundary. 

first,  the  velocity  components  perpendicular  to  the 
boundary  are  set  to  zero  at  the  boundary  and  both  velocity 
components  are  zero  at  the  corners.  The  perpendicular 
velocity  components  just  outside  the  boundary  are  set  equal 
to  the  negative  of  the  perpendicular  components  just  inside 
the  boundary.  The  parallel  velocity  components  just  outside 
t he  boundary  are  set  equal  to  the  parallel  velocity  component 
just  inside  the  boundary.  The  surface  height  iust  outside 
the  boundary  is  set  equal  to  the  surface  height  just  inside 
the  boundary.  The  corners  are  treated  similarly  with  the 
corner  point  being  the  center  of  symmetry  ( e . g .  ,  if  corner 
point  is  (2,2)  ,  the  V(l,2)  -  -V(5,2)  ,  11(1,2)  =  0 

and  U(1  ,1)  =  -  11(3,5)  )  . 

V ej  i  f  i cat  ion 

The  program  based  on  these  difference  equations  and 
boundary  conditions  was  verified  by  duplicating  the  results 
obtained  by  Mader  using  his  SiVAN  code  (Ref  (>)  .  The  problem 
solved  by  Mader  consisted  of  a  .5  meter  radius  hemispherical 
surface  cavity  in  water  three  meters  deep.  The  cavity  was 
the  result  of  a  near  surface  burst.  This  program  used  the 
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computer  simulations  ha  soil  on  t  ho  difference  equations, 
initial  o  ond  1 1  i  on  s ,  anil  boundary  conditions  wore  used  to 
determine  the  sensitivity  of  t  he  tidal  wave  generation. 
Figure  2  illustrates  the  geiieial  shape  of  the  tidal  wave 
generated  by  the  computer  simulations.  The  wave  velocity 
is  33(m/s),  which  is  very  close  to  the  characteristic  wave 
velocity  of  31.3(m/s).  The  initial  condition  used  produced 
a  wave  with  a  period  of  about  SU  seconds  which  is  nearly 
three  times  the  expected  period  of  18  seconds  based  on  the 
relationship  Period  -  11.1  Y*'-1''*  (Ref  3:2"-!).  This 
indicates  the  wavelength  of  l (>()()  meters  produced  by  this 
initial  condition  is  also  about  three  times  too  long.  The 
height  of  the  initial  wave  above  the  surface  was  selected 
as  the  indicator  of  the  wave  generation  sensitivity.  The 
wave  height  is  measured  at  a  nearly  constant  distance  of 
8  kilometers  from  the  explosions.  The  wave  height  was 
obtained  from  the  pr  inter -p lot s  which  expressed  the  wave 
height  within  increments  of  .25  meters.  The  error  bars  on 
the  graphs  indicate  the  variation  in  wave  height  along  the 
plane  tidal  wave.  both  the  linear  and  non-linear  difference 
schemes  used  the  same  time  step  of  0.5  seconds  and  cell 
size  of  50  meters  square. 


Figure  2.  Illustration  of  Plane  Tidal 


Y  t  c  1 d  Sc- ii  sit  i  v  i  I  y 

The  wave  height  sensitivity  to  yield  is  proportional 
to  .  This  does  not  agree  with  the  sealing  of  wave 

height  as  Y^1’  used  for  the  initial  condition  (Kef  5:133). 
However,  it  does  come  close  to  the  scaling  of  wave  height 

r  j  1/4 

as  Y‘J  for  deep  water  and  YJ/  lor  shallow  water  used 
by  Glasstone  (Ref  3:272  273).  The  sensitivity  to  yield  over 
the  range  of  2  -  8  kilotons  is  +4 . 5 (m/ ( kt ) ^ ^ )  . 

The  wave  sensitivity  to  combinations  of  two  yields 
whose  average  is  5  kilotons  is  small.  The  sensitivity  is 
determined  relative  to  the  ratio  of  the  small  yield  over  the 
larger  yield.  As  shown  in  figure  4,  the  wave  height  is  not 
very  sensitive  to  combinations  of  yields  over  a  practical 
range  of  0.5  to  1. 

Spa c  mi g  Sens  i  t  i  v  i  t  y 

The  spacing  sensitivity  was  determined  over  the  range 
of  .6  to  1.4  kilometers.  The  wave  height  did  not  change 
linearly  with  spacing  as  shown  in  I igure  5.  Over  the  range 
of  0.8  -  1.2  kilometers  spacing,  the  approximate  sensitivity 
is  -11.25(m/km)  spacing. 

A  1  ignment  Sens  it  i  v  i  t y 

The  sensitivity  to  m  i  sa  1  ignment  wa*-  termined  for 
misalignments  up  to  .5  kilometer  as  indicated  in  figure  f> . 
The  misalignment  situation  was  set  up  by  having  every  other 
explosion  offset  a  given  distance  from  the  array.  The 
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Figure  3.  Yield  Sensitivity 
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figure  6. 


misalignment  distant  <.  i  .  defined  as  tins  o  I  i  set  distance. 

The  approximate  sen;,  itivitv  over  the  range  ol  0  to  .5 
kilometer  is  -2.7(m/km)  misalignment.  The  misalignments 
did  increase  the  wave  height  variation  along  the  wave  front 
as  indicated  by  the  error  bars. 

I  im  i  in.  Sens  itivitv 

The  wave  height  sensitivity  to  timing  delays  of  0-5 
seconds  is  very  small,  as  shown  in  figure  7.  The  time 
delay  situations  are  modeled  by  having  every  other  explosion 
delayed  the  same  amount  of  t  ime.  In  the  range  of  0  -  2.5 
seconds  delay,  the  sensitivity  is  only  - . 02 h | m/ sec )  time 
delay. 

Bo t  t  pm  S 1 ope  Sens i t  i  i  i t  y 

The  average  slope  of  the  continental  shelf  had  a 
negligible  effect  on  the  wave  height  in  8  kilometers  of 
travel.  The  bottom  was  treated  as  being  smooth  and  with 
constant  slope.  I  he  slope  was  varied  from  0.0  to  0.002 
as  indicated  by  I  igure  S.  The  wave  height  sensitivity  at 
8  kilometers  from  the  explosions  is  +182(m)  per  the  slope. 
Assuming  the  effect  was  constant  over  the  8  kilometers, 
the  sensitivity  was  +  23{in/km)  per  kilometer  traveled  by  the 
wave  per  the  slope.  For  example,  using  the  averaged  conti¬ 
nental  shelf  slope  of  0.0005,  the  expected  increase  in  wave 
height  in  100  kilometers  is  0.7  meters. 


WAVE  HEIGHT  (METERS) 

9.25  9.50  9.75  10.00  10.25  10.50 


Figure  7.  Timing  Sensitivity 
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Non  -Li  near  hi'  fee  t  s 

The  non-linear  effects  are  illustrated  by  comparing 
the  wave  formation  of  a  linear  and  a  non-linear  computer 
simulation.  Wave  contour  printer-plots  of  figures  9  and 
10  show  the  non-linear  and  linear  initial  waves,  respec¬ 
tively,  at  four  kilometers  from  the  explosions.  bach  line 
of  the  contour  plots  indicates  a  change  in  height  of  two 
meters.  The  troughs  or  areas  of  the  water  surface  below 
the  normal  free  surface  are  not  represented  by  contour 
lines,  but  left  blank.  The  plots  show  how  the  non-linear 
effects  produce  a  slightly  smoother  and  larger  wave  in  the 
same  distance.  However,  the  linear  simulation  provides  a 
very  good  approximation  for  the  initial  conditions  of  this 
study . 
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V  .  LOIIC  Ills  I  Oil 


An  array  of  simultaneous  underwater  explosions  is  very 
effective  in  generating  a  plane  tidal  wave.  This  conclusion 
is  based  on  shallow  water  theon  u i t h  the  explosion  represented 
by  a  raised  cone  of  water  which  produces  waves  with  a  wavelength 
of  about  1 6 U 0  meters. 

The  sensitivity  of  the  wave  generation  process  was  deter¬ 
mined  by  the  height  of  the  initial  wave.  The  wave  height  was 
not  very  sensitive  to  misalignments  less  than  200  meters,  time 
del  a vs  less  than  5  seconds,  bottom  slope  less  than  0.0004,  and 
ratios  of  yields  in  the  range  of  0.5  to  1.0.  The  wave  height 
sensitivity  to  spacings  of  1.0  ±  .2  kilometers  is  approximately 
-11.25  meters  per  kilometers  spacing.  The  yield  sensitivity, 
which  depends  on  the  derived  initial  condition,  is  4.55  meters 
per  cube  root  of  the  yield  in  kilotons.  The  non-linear  effects 
increase  the  size  of  the  tidal  wave. 
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V  1  .  Rec omnetula t  i on 

More  research  on  the  initial  condition  is  recommended. 

The  waves  produced  hv  the  raised  cone  of  water  have  a  wave¬ 
length  about  three  times  too  long  and  the  initial  wave  is 
not  preceded  by  a  trough  as  normally  expected  (.Kef  3:275). 

Tn  addition,  the  cone  does  not  conserve  the  volume  of  water. 

A  new  initial  condition  should  be  developed  which  conserves 
the  water  volume  and  produces  waves  with  a  shorter  wavelength. 
It  is  recommended  that  a  cylindrical  cavity  w  i  t  It  its  depth 
equal  to  its  radius  (l no  meters)  and  with  a  lip  of  water 
surrounding  it  be  used  as  a  new  initial  condition. 
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This  report  considers  an  array  of  simultaneous  underwater 
nuclear  explosions  generating  plane  tidal  waves  on  the  continental 
shelf.  Computer  simulations  are  used  to  study  this  wave  genera¬ 
tion  process.  They  are  based  on  non-linear  shallow  water  theory 
with  an  initial  condition  of  a  stationary  raised  cone  of  water. 

The  sensitivity  of  the  wave  height  to  different  yields,  combina¬ 
tion  of  different  yields,  spacing,  alignment,  tuning,  and  bottom  ^ 
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slope  is  determined.  The  reference  values  used  (or  this  study 
are  a  yield  of  five  kilotons  INI  ,  one  kilometer  spacing,  water 
depth  of  100  meters,  and  no  bottom  slope.  In  addition  t o  the 
sensitivities,  the  amount  of  non-linear  effect  is  illustrated 
by  comparison  of  wave  contour  p r  i  n t e r  -  p 1 ot s  generated  by  linear 
and  non-linear  shallow  water  theory.  The  computer  results 
indicate  the  wave  generation  is  relatively  insensitive  to  real¬ 
istic  variations  in  the  parameters. 
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